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ABSTRACT 

Recently, HESS and other air Cerenkov telescopes have detected a source of TeV 7-rays coincident with 
the Galactic center. It is not yet clear whether the 7-rays are produced via leptonic or hadronic processes, 
so it is important to consider possible acceleration sites for the charged particles which produce the 7-rays. 
One exciting possibility for the origin of these particles is the central black hole, Sgr A*, where the turbulent 
magnetic fields close to the event horizon can accelerate protons to TeV energies. Using a realistic model of the 
density distribution in a 6 pc x 6 pc x 6 pc cube at the Galactic center, we here calculate the trajectories followed 
by these TeV protons as they gyrate through the turbulent medium surrounding Sgr A*. Diffusing out from the 
black hole, the protons produce TeV 7-rays via ir° decay following a collision with a proton in the surrounding 
medium. After following over 222,000 such trajectories, we find that the circumnuclear ring around Sgr A* 
can reproduce the observed 0.1-100 TeV HESS spectrum and flux if the protons are injected into this medium 
with an effective power-law index w 0.75, significantly harder than the observed photon index of 2.25. The 
total energy in the steady-state 1-40 TeV proton population surrounding Sgr A* is inferred to be w 5 x 10 45 
ergs. Only 31% of the emitted 1-100 TeV protons encounter the circumnuclear torus, leaving a large flux of 
protons that diffuse outward to contribute to the Galactic ridge emission observed by HESS on scales of > 1°. 
Subject headings: acceleration of particles — Galaxy: center — radiation mechanisms: nonthermal — gamma 
rays: theory 



1. INTRODUCTION 

The Galactic center is a complex environment, harboring 
what is believed to be a supermassive black hole, Sagittar- 
ius A* (Sgr A*), surrounded by a blend of new and old 
stellar populations, thermal and nonthermal gas components, 
and a wide array of both compac t a nd diffuse hi gh-energy 
sources (see iMelia & Falckal2001l or IMelia 2007 for recent 
reviews). In an exciting recent development, the Galactic cen- 
ter has been identified as a source of TeV 7-rays by the air 
Cerenkov telescopes HESS dAharonian et all 12004 120 06b), 
Whip ple dKosack et alJl2004. CANGA ROO (iTsuchiya et all 
120041) and MAGIC (lAlbert et al.ll2006l) . While many of the 
original observations initially measured diff erent spectra and 
lumin osities for the Galactic center source dAharonian et al.1 
2004), subsequent reprocessi ng and analysis has brought all 
the detections into agreement (Katagiri et al. 2005). Here, we 
concentrate exclusively on the HESS data as they provide the 
best constraints on the source properties. 

The signal from the Galactic center was detected by HESS 
in observations conducted over tw o epochs (June-August 
2003 a nd March-September 2004; lAharonian etail 12004 
2006b:) with a combined ~ 3 8<r excess above the back- 
ground (lAharonian et al.ll2006bl) . The spectrum of this source 
is a power law with photon index 2.25 ± 0.10, and the to- 
tal flux above 1 TeV i s (1.87 ± 0.3 0) x 10 8 m" 2 s" 1 
( lAharonian et al.ll2006bT) . As argued bv lCrocker et all d2005l) . 
this HESS source does not appear to be coincident with 
the EGRET Galactic center source 3EG J 1746-2851. There 
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are two reasons for this: first, the EGRET source ex- 
cludes the Galactic center at the 99.9% confidence level 
dHooper & Dingusll200H IPohl 120051) . In contrast, the HESS 
source is coincident within ~ 30" of Sgr A*, though with 
a centroid displaced roughly 7" (~ 0.4 pc) to the East of 
the Galactic center (Aharoni an et al.l l2006bl) . Second, the 
EGRET spectrum extrapolated into the HESS energy range 
over-predicts (by a factor ~ 20) the TeV 7-ray flux of the 
Galactic center source. Although it is possible that both 
sources are associated with the SNR Sagittarius A East, the 
HESS TeV source is probably associated with the black hole 
itself, given its angular proximity to Sgr A*. 

If the TeV 7-rays do originate from Sgr A*, then this 
presents somewhat of a theoretical puzzle: given what 
we understand about Sgr A* and the interactions with its 
nearby environment, this object should not be directly pro- 
ducing a significant flux of TeV photons. Several pos- 
sible scenar ios for producing these 7-rays have been pro- 
posed (e. g.. lAtovan & Dermer 2004; Aharo nian & Neronovl 
2005altt [Hinton & Aharonian! 120061) . including that of 
Liu et al. J2006I) . These last authors show that protons can 



be energized to TeV energies by stochastic acceleration in 
a magnetically dominated funnel close to the black hole (~ 
2 — 30 Schwarzschil d radii). The proton acceleration models 
of lBeckeret"ai] (2006) may also be active (or even dominant) 
in Sgr A*. In either case, the relativistic protons will then 
diffuse along magnetic fields from small scales to large radii 
without experiencing adiabatic loss as is the case for ejections 
of plasma blobs. As the protons diffuse out through the sur- 
rounding medium, they may scatter with hydrogen nuclei in a 
shocked stellar wind region, and with the circumnuclear disk 
surrounding the black hole. Such pp scattering events produce 
7r°'s, which subsequently decay into two photons. Thus, this 
model suggests that protons random-walking their way into 
the gaseous few central parsecs surrounding Sgr A* may ulti- 
mately be the source of the H ESS -detected TeV 7-rays (e.g., 
lAharonian & Nerono v 2005b). It is this process that we aim 
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to examine carefully in this Letter by carrying out a highly- 
detailed analysis of the proton transport and interaction, start- 
ing at the acceleration site near the black hole and extending 
out as far as several parsecs. We assume a distance to the 
Galactic center of R = 7.94 kpc (lEisenhauer et alJl2003l) . 

2. CALCULATIONS 

2.1. Outward Proton Diffusion from Sgr A * 

The motion of a relativistic proton moving through the ISM 
is governed solely by the Lorentz force due to the interstellar 
magnetic field: dv/dt = (v X fl) /7, where 7 is the Lorentz 
factor for a proton with velocity v and S~2 = eB /mc is the 
gyrofrequency for a particle with charge e and mass m under 
the influence of a magnetic field B. Thus, given a descrip- 
tion for the magnetic field in the computational volume, the 
trajectory of the proton can be followed exactly. As there is 
no net electric field, the magnitude of the velocity vq remains 
constant. 

Proton trajectories are calculated in a 6 pc x 6 pc x 6 pc cube 
consisting of 10 6 equally spaced cells centered on the Galactic 
center. The hydrogen density nn within this volume i s take n 
from the results of simulations by iRockefeller et alj d2004l) . 
These authors computed the density distribution in this region 
caused by the interactions of stellar winds from the surround- 
ing young stars. In addition to the stellar wind gas, the vol- 
ume also contains a high-density 'torus' with an inner radius 
of 1.2 pc and a thickness of 1 pc representing the observed 
circumnuclear disk containing molecular gas. 

The magnetic field is assumed to be generated everywhere 
in this region with an intensity that is proportional to njj ■ The 
average density in the model stellar-wind gas (taken to be any 
region where jih < 3 x 10 3 cm~ 3 ) is {n™) = 121 cm~ 3 , 
while it is (rtg'} = 233,222 cm~ 3 within the model molecu- 
lar torus. Taking kT =1.3 keV as the aver a ge temperature of 
the st ellar-wind gas (B aganoff et alj|2003t : IRockefeller et alj 
120041). 100 K for the temperature of the molecular torus 
(IRockefeller et al J 120 04 and references therein), and assum- 
ing equipartition, the average field intensity Bq is w 3 mG in 
the stellar-wind region and « 0.35 mG within the torus. As 
the proton moves through different densities in the computa- 
tional domain, Bq is scaled to the appropriate value, indicated 
by the local gas density. 

While the above procedure scales the magnetic field in- 
tensity at any position in the grid, we also require a physi- 
cal description of the field direction at each point. A use- 
ful prescription f or describing a turbulent magnetic field 
was developed by Giaca lone & JokiplH (1 19941) to analyze the 
fundamental physics of i onic motion in the ISM (see also 
iKowalenk o & Melia 2000). This quasi-analytic approach as- 
sumes that the magnetic field fluctuations are static and follow 
a Kolmogorov spectrum (though in principle any type of tur- 
bulence may be modeled in this way). At every position in 
the grid, the three-dimensional magnetic field is written as the 
sum over a given number of transverse waves of random po- 
larization, with a wavevector oriented in a random direction. 
In contrast to iGiacalone & JokipiH d 19941) . we add together 
fluctuations from 200 wavelengths, ranging from O.l^o/^o 
to lOvo/flo, where Slo is the proton gyrofrequency at Bo. In 
this fashion, we may then generate all three components of 
the magnetic field as functions of the position coordinates, for 
each proton released from the origin with initial velocity vo 
in a random direction. 

With the above ingredients, the trajectory of each proton 



can be calculated exactly as it traverses the volume around Sgr 
A*. The calculation is completed once the proton leaves the 
computational domain. Since the observed-frame time step is 
proportional to 7 2 , computational resources limited the cal- 
culations to proton energies between 1 and 100 TeV. A total 
of 222,617 proton trajectories were calculated with energies 
uniformly distributed in this range. Figure Q] shows a typical 
trajectory calculated in this fashion. 

2.2. The Pion Production Rate and Gamma-ray Emissivity 

As the proton random-walks its way toward the edge of the 
calculation domain by scattering off the turbulent magnetic 
field, it may collide with a low-energy proton in the ambient 
medium and produce pions via the reaction pp — > ppir°Tr , 
where a variable number of pions are produced to conserve 
charge. The neutral pion will subsequently decay into two 
photons. In principle, one or more of the two protons in each 
scattering event may retain (or gain) sufficient energy to pro- 
duce additional pions in subsequent interactions, but we will 
ignore their contribution to the overall pion production rate 
for this application. In addition, a charge exchange may occur 
(roughly 1 /4 of the time), in which a neutron is produced in 
the exit channel. A fit ting formula for the cro ss-section of pp- 
scattering is given by Eidelma n et alj (|2004), and is typically 
a vv ~ 40 mbarns for proton energies between 1 and 10 TeV. 

To compute the fraction of relativistic protons that undergo 
pp-scattering, we first split the trajectories into 21 energy bins 
(log(-EyeV) = 12, 12.1,..., 13.9, 14). Each proton path is 
then followed and if the particle enters a volume element with 
tt-h > 3 x 10 3 cm~ 3 (these regions will dominate the 7-ray 
emissivity) the incremental optical depth dr — n-na pp (E)dl 
to pp-scattering within this grid cell is calculated. Here, dl is 
the distance the proton traveled since the previous time step 
(w 1 /Vl) and fin is the average density encountered by the 
proton as it moved a distance dl. As long as the particle 
remains in that cell, we continue to sum dr. If the proton 
scatters out of the volume element and then subsequently re- 
turns, the new values of dr are added to the previous total. 
This procedure is followed for each proton trajectory in each 
of the 21 energy bins. The fraction of protons at energy E 
that undergo a pp-scattering event in a volume element with 
n H > 3 x 10 3 cm~ 3 is f(E) = 1 - e r^ dT ^ N , where N is 
the number of calculated trajectories at energy bin E. 

A power-law spectrum of proton energies in the relativis- 
tic regime is a natural cons equence of many acceleration pro- 
cesses (e.g jLiu et al.l2006l) . but the exact value of the spectral 
index depends on the details of how the protons escape. As we 
shall see below, the rate of diffusion and subsequent scatter- 
ing of the protons is highly energy dependent. Therefore, we 
treat the "injected" spectral index as a variable to be fixed by 
the fitting procedure. By "injected" we mean the distribution 
of protons leaving Sgr A*'s region of influence and entering 
the wind-shocked region surrounding it. 

The proton spectrum is thus written as dn/dE^ = 
K(E/E m i n )~ a , where K oc £, the total energy in protons 
in ergs between E ni i n = 1 TeV and 40 TeV. The spec- 
trum of scattered protons in a grid cell is thus dn/dE = 
(dn/dEi n j) f(E), where f(E), defined above, is the fraction 
of protons with energy E which collide with an ambient pro- 
ton in that volume element. However, f(E) is only known 
down to log-E = 12, so dn/dE is extrapolated to 0.1 TeV 
using the chosen power-law index a and the value of f{E) 
at \ogE — 12. As 7-rays at a given energy are preferen- 
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FIG. 1. — (Left) Complete trajectory of a log(i?/eV) = 12.4 proton from ejection at Sgr A* (the origin) until it escapes the computational volume at 2 = 3pc. 
The x-y plane is the sky plane. (Right) The magnetic field strength and gas density encountered by this proton on its trajectory. The particle traveled for 1,211 yrs 
and traversed 371.3 pc on its journey outward from Sgr A*. It spent most of its time in the lower-density stellar wind gas but did have numerous encounters with 
the high density molecular torus. The dotted line in the upper-panel indicates the lowest magnetic field strength for which the Lorentz equation was solved. For 
the extremely low density zones which had magnetic field strengths below this limit the particle was moved in a straight line until it encountered higher density 
gas. 

rived value of a — 0.75 is much harder than the "standard" 
power law with index ^2 — 2.5 that one might expect from 
simple first-order Fermi scattering. Perhaps this may be taken 
as an indication that Sgr A* is not the source of energetic 
hadrons producing the 7-ray glow in its vicinity. Alterna- 
tively, the intrinsic a ~ 2 proton spectrum may have been 
altered by interactions much closer to the black hole. Another 
possibility is that the required hard proton spectrum might be 
an indication that stochastic acceleration continues as the par- 
ticles diffuse outwards from the initial injection region. Fu- 
ture work will investigate the origin and viability of the hard 
proton spectrum. 

We find that only about 30% of the computed trajecto- 
ries encounter the circumnuclear torus; the rest travel unim- 
peded to larger distances. However, this fraction is highly 
energy dependent: ~ 73% of the emitted protons encounter 
the torus at log(E/eV) =12-12.4, but this drops to 47% at 
log(J5/eV) = 13, and finally to 5% at log(£7/eV) = 14. As 
a result, the predicted 7-ray spectrum exhibits a roll-over at 
energies > 20 TeV (Fig. g). This is a robust prediction of 
models where Sgr A* is the source of the relativistic protons 
and can be tested with future HESS observations. 

The energy-dependent proton interaction rate also suggests 
that a significant flux of TeV protons escapes from the in- 
ner few pc surrounding Sgr A* to interact with more distant 
molecular gas. This would be consistent wi th the observed 
emissio n from the Galactic center ridge. Ah aronian et alJ 
(2006a) report that the inner 200 pc of the Galaxy glows 
in TeV 7-rays, with the dominant emission correlated with 
the distribution of dense molecular gas. In addition, the 
spectrum of the Galactic ridge emission is observed to be 
close to that of the Galactic center source. The best-fitting 
model for the Galactic center HESS source predicts an en- 
ergy of £ = 5 x 10 45 er gs in protons between 1 and 40 TeV. 
Aharo riian et alJ d200 6a) estimate that a hadron distribution 
with an energy ~ 10 49 ergs (from 4^40 TeV) is required to 
produce the observed Galactic ridge emission. The ridge re- 
gion is much larger than the volume within which the Galac- 



tially produced by protons with ~ 10 x greater energy, this 
extrapolation will not greatly affect the predictions between 
0. 1 and 1 TeV. The fraction of protons that interact with the 
circumnuclear torus becomes very small at energies greater 
than 100 TeV (see § [3j, so it is not necessary to extrapolate 
dn/dE beyond 100 TeV. Now armed with the dn/dE profile 
within each cell, we calculate the tt° and 7-ray emissivity for 
different values of £ and a using the formalism described by 
ICrocker et alJ ((2005). The emissivities are then converted into 
fluxes before the values are summed to produce the predicted 
7-ray spectrum from 0.1-100 TeV. 

3. RESULTS 

We compared the model spectra to the observed 2003 
and 2004 HESS data between 0.3 and 30 TeV (taken from 
lAharonian et all 12004 12006b) using x 2 fitting. The low- 
est x 2 was found for £ = 5 x 10 45 ergs and a = 0.75 
[x 2 /d.o.f= 41.7/31 (29.1/31) using the smaller (larger) of 
the observed error bars, and d.o.f. = degrees of freedom] . This 
value for the spectral index is very robust with A\ 2 ~ +5 for 
a = 0.7 or 0.8. The predicted 0.1-100 TeV spectrum and 
1-10 TeV image of the Galactic center with these parame- 
ters are shown in Fig. [2] The expected 1-10 TeV luminosity 
from this spectrum is 4.7 x 10 34 erg s _1 . The predicted 1- 
100 TeV flux is 1.83 x 10~ 8 m~ 2 s _1 , in good agreement 
with the observed value of (1.87 ± 0.30) x 10~ 8 m~ 2 s _1 
( lAharonian et alj|2006bl) . 

The fact that the best fit corresponds to a = 0.75 indicates 
that the spectrum of injected particles must be significantly 
harder than the observed TeV 7-ray spectrum (for which the 
spectral index is 2.25 ± 0.10). This is simply because the 
higher energy protons diffuse outward from the Galactic cen- 
ter faster than the low energy ones and as a result traverse, 
on average, a much smaller path length. These particles are 
then less likely to encounter the molecular torus and undergo 
a pp-scattering event that will produce a TeV 7-ray. A harder 
intrinsic spectrum is therefore required to compensate for this 
deficit and account for the observed spectral shape. The de- 
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FIG. 2. — (Top) The solid line is the predicted spectrum from the 
Galactic center region due to pp-scattering. The total energy required in 
protons between 1 and 40 TeV is 5 X 10 45 ergs, and the injected spec- 
trum is dn/dEinj oc E~ °- 75 . T he open data points ar e the 2003 ob- 
served values from HESS taken from Aharonian et al. 12004), while the solid 
points are from the 2004 HESS observations (Aharonian et al. 2006b). (Bot- 
tom) Contours of 1-10 TeV flux predicted from the same model that pro- 
duced the spectrum in the left-hand panel. The contours are at log / = 
-15, -14, -13, -12, -11, -10, and -9 [s^ 1 m~ 2 ]. The oblateness evi- 
dent in this image is due to the tilt angle of the molecular torus. Sgr A* is at 
the origin. 

tic center TeV source is produced, so the escaping protons 
would build up the much larger energy content over time. The 
diffusion time to reach such large distances is > 10 4 years 
dAharonian et al.l r2006al) . so the particle accelerator would 
have to be at least that old. 



Our simulations show that an injected proton index a w 
0.75 is required to produce the observed 7-ray index a bs ~ 
2.25, a steepening of 3/2. Homogeneous diffusion with a 
Kolmogorov or Bohm turbulence spectrum would produce a 
steepening of 1/3 and 1, respectively. Our treatment in this 
paper has followed the diffusion of protons taking into ac- 
count the variable 3-D gas density and magnetic field sur- 
rounding Sgr A*. On the other hand, diffusion calculations 
assume a homogeneous medium, which produces misleading 
results when applied to regions as complex as that near Sgr 
A*. As such, our finding that a steepening of 3/2 is required, 
instead of ~ 1, may be traceable to our more accurate treat- 
ment of inhomogeneity. 



4. CONCLUSIONS 

The HESS observations suggest that there may be a com- 
mon cause for both the Galactic-center point source and the 
diffuse emission along the ridge. TeV protons accelerated 
near Sgr A* can account for the observed flux and spectrum 
of the central source through scattering events with ambient 
protons within the circumnuclear ring. The proton spectrum 
injected into the gas surrounding Sgr A* must have a hard 
power-law index of 0.75, implying a non-standard accelera- 
tion process or interactions close to the acceleration region. 
Our simulations have resulted in a predicted 1-10 TeV im- 
age of the source for comparison with future high-resolution 
observations. We have found, in particular, that the predicted 
TeV image mirrors closely the structure of the torus, including 
an oblateness due to projection effects in the plane of the sky. 
The total energy required in protons between 1 and 40 TeV 
is 5 x 10 45 ergs and a significant fraction of the protons do 
escape the neighborhood of Sgr A* without undergoing a pp 
collision; these presumably diffuse to much larger distances 
where they can interact with other molecular material, possi- 
bly accounting for the observed Galactic ridge emission. Al- 
though future calculations are required to confirm this sug- 
gestion, and o ther origins for the observed HESS source ar e 
possible (e.g., IWang et al J [20061: iHinton & Aharonianll2006l) . 
this work supports the view that Sgr A* may be an important 
site for particle acceleration at the center of the Galaxy. 
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